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O 

Cn ' In this work we study, on a finite and periodic lattice, a class of one-dimensional (bimolecular 

O „ and single-species) reaction-diffusion models which cannot be mapped onto free-fermion models. We 

Q^ ■ extend the conventional empty-interval method, also called interparticle distribution function (IPDF) 

^0 ' method, by introducing a string function, which is simply related to relevant physical quantities. 

As an illustration, we specifically consider a model which cannot be solved directly by the conven- 
»vj ' tional IPDF method and which can be viewed as a generalization of the voter model and/or as an 

epidemic model. We also consider the reversible diffusion-coagulation model with input of particles 
and determine other reaction-diffusion models which can be mapped onto the latter via suitable 
1-^ ' similarity transformations. Finally we study the problem of the propagation of a wave-front from 

an inhomogeneous initial configuration and note that the mean-field scenario predicted by Fisher's 
equation is not valid for the one-dimensional (microscopic) models under consideration. 
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C^ . I. INTRODUCTION 

fi: . . . 

' . ' Reaction-diffusion models (RD) play an important role in the description of classical interacting many-particle non- 

(-H ' equilibrium systems and have been extensively investigated in the last decade KM- Often these systems have been 

Q treated by mean-field techniques (e.g. rate equations) which give rise to non-linear partial differential equations (such 

O [ as, e.g., the Fisher equation ||]). The latter represent difficult mathematical problems: e.g. the Fisher equation cannot, 

in general, be solved exactly. The mean-field methods can accurately describe the behaviour of RD systems in higher 

dimensions, where the correlations do not dramatically change the physics of the models. However, in one spatial 

^ , dimension where the ffuctuations play a crucial role, these mean-field treatments fail. In this sense, a satisfactory 

f^^ ' understanding of RD models in lower dimensions requires exact solutions, which are scarce. In some cases, however, 

C^l ! some RD models are known to be solvable. These cases can essentially be classified into four categories: (i) models 

[^ • for which the equations of motion of correlation functions are closed M ; (ii) the free-fermion models M (or systems 

which can be mapped onto the latter, see H,^); some other one-dimensional RD models can be solved by the (iii) 

Matrix Ansatz method B first introduced to study the steady-states of the asymmetric exclusion process (ASEP) 

r-^ ' and which has been extended to other multispecies RD models where the total number of particles is conserved. A 

"••-y . dynamical version of the matrix Ansatz [Q has also been proposed to study the dynamical properties of the models 

C^ ■ for which the equations are closed (on periodic as well as open chains). Some other one-dimensional models can be 

C , solved by (iv) the empty- interval method, also called interparticle distribution function (IPDF) method pd|-p^ , first 

I ■ introduced for the study of the diffusion-coagulation model. The solution of various one-dimensional RD models have 

' O ] been obtained from the diffusion-coagulation model via similarity transformations [p|pq] . It has been established that 

^ . the latter solvable models correspond to free-fermion systems [0| . 

The purpose of this work is to present a generalization of the IPDF method and to apply this technique to solve a 
class of one-dimensional stochastic models which cannot be mapped onto free-fermion systems. In fact much attention 
. !^ ^^^ been given to free-fermion systems in various contexts (using fermionic algebra ||7| or via the traditional IPDF 
k> ' method, in the continuum limit |iy,f2| as well as on discrete lattice [nsllSl). The situation is different for the models 
; I ' considered here, for which only a few exact results are known. 

Cd , The paper is organized as follows: in the next section we briefiy recall the formalism which we employ. In section 

III we introduce the string-function which is the key to our analysis and determine the constraints necessary to have 
solvable situations. In section IV we solve the general equations of motion for the string-functions of reaction-diffusion 
models which cannot be mapped onto free-fermion systems. The latter provides the exact expression of the density 
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and the instantaneous nearest-neighbour (two-point) correlation function. We also present an approximative, and, 
recursive scheme to compute the (other) instantaneous two-point correlation functions. In section V, we specifically 
consider a model with branching and coagulation reactions which cannot be solved directly by the traditional IPDF 
method. In section VI, we solve a reversible diffusion-coagulation model with external input of particles. In section VII, 
we take advantage of the results of the previous section to solve other related models via similarity transformations. 
In section VIII we study, for the models introduced and solved in section V and VI, the problem of the propagation 
of a wave-front starting from an inhomogeneous initial state and observe that the mean-field scenario predicted by 
Fisher's equation ^j^j^^ is not valid at the microscopic level. Finally, the section IX is devoted to the conclusion. 

II. THE FORMALISM 

Before generalizing the IPDF method, it is useful to briefly review the so-called stochastic Hamiltonian formalism . 

It is known that models of stochastic hard-core particles are soluble on some manifold on wich the equations of 
motion of their correlation functions close [^ . In this work, we concentrate on one-dimensional bimolecular single- 
species reaction-diffusion systems. 

Consider a periodic chain with L sites (labelled from 1 to L). On the lattice, local bimolecular reactions between 
single-species particle A, with hard-core, take place. Each site can be empty (denoted by the symbol 0) or occupied 
at most by a particle of type A denoted in the following by the index 1 . The reactions occuring on the sites j and 
j -|- 1 are specified by the transition rates, which here are assumed to be site-dependent^ according to F^^, where 

a, 13,-1,5 = 0, 1: \/{a, /3) 7^ (7, (5), F^^ : a -f /3 — > 7 + (5. Probability conservation implies F^^ = - Y.(a.p)^(a' ,0') ^Ip 

andF2;'^>0, V(a,/3)^(7,<S) 

For example the rate F}? corresponds to the reaction AA — > A^ and F}} = -(F}? + F?} + F?J). 

The state of the system is represented by the ket \P{t)) = X)/n> ^(i"-}'^)!^*)' where the sum runs over the 2^ 
configurations. At site i the local state is specified by the ket \ni) = (1 0)"^ if the site i is empty, \ni) = (0 1)^ if the 
site i is occupied by a particle of type ^ (1) . 

It is by now well established that a master equation can be rewritten formally as an imaginary time Schrodinger 
equation: ■^\P{t)) = —H\P{t)), where H is the Stochastic Hamiltonian which governs the dynamics of the 
system. In general, it is neither hermitian nor normal. Its construction from the master equation is a stan- 
dard procedure (see e.g. |l|,^) The evolution operator H = X]i=i-^jj+i ^^ts locally on two adjacent sites, with 
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— iJjj+i = I ^10 y^io pio pio I : where the same notations as in reference [j4|-|6|] have been used. Probability 

conservation implies that each column in the above representation sums up to zero. 

The left vacuum (x|, which is defined as (x| = Sinld*^}!' locally has the representation (x| = (1 1) (8) (1 1) with 
the property: (x|i?jj+i = 0. 

Below we shall assume an initial state \P{0)) and investigate the expectation value of an operator O (observables 
such as density etc.): {0){t) = (x|Oe~^*|P(0)). For general s— species bimolecular reaction-diffusion systems, there 
are (s-|-l)^ possible rates which have to fulfill the (s-j-l)^ probability conservation constraints. Thus general s— species 
bimolecular reaction-diffusion systems are characterized by (s -|- 1)^ — (s + 1)'^ independent parameters P,p[. If one 
imposes to these parameters 2s^ appropriate constraints, the equation of motion of correlation functions close and 
the system is formally soluble in arbitrary dimensions. Here, we focus on the case s = 1, and so we have 16 — 4 = 12 
independent rates and 2 closure constraints. 

For single-species bimolecular processes, with the notations: 

pOO , pl0\ 

-■- 11 + J- 11; : 

r?; + F?}), (1) 

the closure constraints are the following [H: 

D, = i}^T\^ + Til (r?? + r?}) = pi? + Tl\ (F?8 + T%) , 

i?i = ^ F8J + Til - (r?? + r}?) ^ Til + r\l - (r™ -h fJ?) , (2) 

With help of the relationships M : 
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- {rimHrn-l^rn) = M + Bi{nra-l) ~ Ci (n„) +Di(n™_inm) , 

- {nmHra.m+l) = ^2 + -B2(nm+l) - C'2(n,„) + i:)2 ('T-m"-m+l) (3) 

The equation of motion of the density at site m (on a periodic chain) reads: 

-rMm){t) = -(xl«me-^*|P(0)) =^1+^2+ Bi(n„_i)(i) + S2(n„+i)(i) 
at at 

- (Ci + C2)(n„)(i) + L>i(n,„_in„)(i) + i:>2(rimnm+i)(i) (4) 

In order to iUustrate the physical meaning of the models studied in this work, let us consider the latter at the 
mean- field level (in the continuum limit), i.e. we assume first: {uxit)) -^ pMF{x,t) and {nxnx±i){t) — {pMF{x,t))^. 
At this level of approximation, we note (see the section VIII) that the equation of motion (H) of some microscopic 
reaction-diffusion models studied in this work, provided that Di < and D2 < 0, is a non-linear partial-differential 
equation of Fisher-type P Jl^ , p^ : 

d ^ d"^ ^ _ _ 

jr;PMF{x,t) = B—-^pMF{x,t) + kipMF{x,t) - k2{pMF{x,t)f , (5) 

with pMF{x,t) = pMF{x,t) - (f). When Ai = A2 = A ,Bi = B2 = B ,Ci = C2 = C, and Di ^ D2 = D < 0, we 
have 2D(j> = {C - B) + ^{B-Cf -AAD, fci = 2^{B-CY -4,AD > which denotes the constant describing the 
growth and k2 = —2D > is the constant describing the saturation according to the local dynamics ||3|, |l6| , [l7| . 

Fisher equation (0) admits two homogeneous steady-states, namely: pmf{x, 00) = pmf{oo) — ki/k2, which is 
stable and another, unstable, steady-state: pj\/f(oo) = 0. 

Although Fisher equation (pSJ) cannot be solved exactly, it is known |1Q,0 ^^^^ ^^^ approach towards the steady- 



state from inhomogeneous initial states (e.g. pmf{x,0) — |^8(xo — x), where Q{x') = 1 if a;' > and Q{x') — 
otherwise) is characterized by a wave-front: pMF{x,t) ~ f{x — ct) propagating with a celerity c > 2^/kiB and 
satisfying the non-hnear differential equation: Bj^f{z) -\- Cj^f{z) -\- kif{z) — k2{f{z)Y = ]l6| |. 

In this work we obtain exact expression for the density from the A'^— body description of some reaction-diffusion 
models (on finite and periodic lattice), for which Di = D2 ^ and for which, the mean- field description in the 
continuum limit is given by non- linear partial-differential equation of Fisher-type (0). Therefore, with help of the 
(microcopic) exact results obtained in section V and VI, we are able, in section VIII, to discuss the validity of the 
mean-field description a la Fisher. 

III. THE STRING FUNCTION 

In this section we introduce the quantity which is the key to our analysis, i.e. the string function Sx^y{t). We also 
derive the constraints for which the equation of motion of the latter is a closed hierarchy. In the sequel, we solve the 
latter providing the density of particles and some two-point correlation functions. 



Instead of considering the standard "empty- interval function" |11-15|, we focus here on the more general form 
(1 <a; < y<x + L): 

Sx,y{t) = {{a - hnx)(a- hnx+i) . . . (a - bny^2){a- bny^i)){t), (6) 

where a and b are non- vanishing numbers. 

This expression reduces to the empty-interval function when a = 6 = 1 |0-E5l. Hereafter we will derive the 
equation of motion of the quantity Sx,y{t) and determine which constraints are necessary and sufficient to close the 
latter. An alternative approach would consist in considering the empty-interval function (with a = 6 = 1 in (jg)) 
and obtaining a solution of some related reaction-diffusion model via a similarity transformation. This approach has 
been extensively studied for the free-fermion models [ [l3|Jl5[| where solutions of RD systems are obtained from the 
solution of the (free-fermion) diffusion-coagulation model. The latter approach is investigated in section VII. For 
l<x<y<x-\-L, the equation of motion of Sx,y{t) reads: 

-7:Sx,y{t) = - ((a - bnx)Hx-i^x{a - bux+i) ... (a - buy^i)) (t) - {{a - bnx){a - bux+i) ... (a - bny^i)Hy^i^y) (t) 
y-2 
- 51 ^^" ^ ^"^)(" ^ bux+i) ... (a - bny^i)H.j,j+i) (t) (7) 

j=x 



If the following five constraints are fulfilled, the dynamics of Sx.y{t) is described by a closed hierarchy of equations: 



(1) 
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The interesting point is that this five constraints are generally independent of the previous ones (0). Therefore, in 
general, models which are solvable via the approach presented here are not on the 10— parameters manifold described 
by (0) where the equation of motion of correlation functions are closed. 

When the constraints (^ are fulfilled, the equation of motion of the string-function Sx,y{t) on a periodic lattice of 
L sites is the following (for aCi ^ bAi and aC2 ^ 6^2) n: 



f ffSx.,y{t) = (aCi - bAi)Sx+i,y{t) + (aC2 - bA2)Sx,y-i{t) - ^5.-i,,(i) 



{B,+B2+C^+C2)Sx,y{t)+ Tll + {^){Tll+Tll) + {^) Til {v ^ x - l)Sx.y{t); [l < x < y < x + L) 
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where the boundary condition Sx,x{t) = 1 is obtained from the requirement that 4:Sx,x+i{t) = ~i>4:{nx{t)) ■ 

Let us note that when a = 5 = 1, we recover the usual constraints of the IPDF method: V^l = VW, Fqi = r}i, 

pOO _ pOO _ n 
J^ 01 ^ ^ 11 ~ ^■ 
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IV. SOLUTION OF THE EQUATION OF MOTION OF THE STRING-FUNCTION 

Equations of motion for the string-function have been intensively studied for frce-fermion systems both in the 
continuum limit ||ll|,|l^ and on discrete lattices (periodic and open boundary conditions) ||l^-[l5t]. However, for systems 
which cannot be mapped onto free-fermion systems, a few results have been obtained: Doering and ben-Avraham 11^] 
have obtained the stationary concentration and the relaxation spectrum of a reversible diffusion-coagulation model 
on an infinite lattice in the continuum limit. Later, Peschel et al. [M have studied the relaxation spectrum of the 
Fourier transform of the string-function on a periodic lattice with help of the conventional IPDF method. 

It is useful to introduce the following notations: 
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^If aCi = Ml and aC2 = M2 
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and tlie equation of motion of Sx,y reads: 



x,y-\-l 



(t) - (7 + S)Sx,y{t) -{y-x- l)SS:,,y, {l<x<y<x + L) 



■J^Sx,x + L{t) = —L5Sx,x + L{t) 



,where 7 -I- 5 = Bi -I- B2 -I- Ci -I- C2, and S = - 
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An example of model which dy- 



namics is described by such a system of equations (with a = 6) is the random sequential adsorhtion (RSA) process of dimers: 



Herefater we solve the equation of motion (g) for the case S ^ (S — corresponds to the free-fermion case), Di ^ 
and I?2 7^ Oq with the additional condition, a2 — ai = a ^ and /3i = /32 = /? ^ 0, which corresponds to the 
restriction to unbiased systems. 

We also introduce the following auxiliary function: TZx.y{t) = ii^'^Sx.yit), where we choose fi = ^/^ = i/§ and 

q — ^/aiP2 = \fo^ and, we solve (P) for the general case where g y^ (Note that /3i = ^^ 1?^ = i?i = g = 0, with 
i € (1, 2); see footnote (2) ). With these notations, the equation of motion (g) becomes: 

ft'^x,y{t) = § Ee=±i {^^+e,y W + T^x,y+e{t)} - inx,y(t) - (y - a;)(57^,,^, (t) ; {l<x<y<x + L) 
itTl,,x+L{t) = -LSn^.x+Lit) (11) 

7^x,x(i) = 1, 

The stationary solution of the system ([lj)is obtained with help of the properties of Bessel functions of first and 
second kind, respectively, J^{z) and Yy (z) M. In fact, the structure of the equation ( pi] ) for 1 < a; < y < L suggests 
the Ansatz: Rx^y{oo) = A^Jy^^+i^i'^q/S) + B ^Yy^x+uii'^q / 5) . Inserting this expression into (pl|), we obtain uj = 7/(5. 
Therefore we have: 

i?,,j,(oo) = ALJy_,+^/s{2q/5) + BLYy_.,+^/s{2q/S) (12) 

The quantities Al and El are determined with help of the boundary conditions: Rx,x{oo) = 1 = ALJi+-y/s{2q/S) + 
BLYi+^/s{2q/5) and Rx^x+l{oo) ^ ALJL+-,/s{2q/5) + BLYL+^,s{2q/5) = 0. It follows 

which, with (O), provides the stationary expression of TZx.y{oo) and thus we obtain the stationary expression of the 
string- function: 



Sx,y (OO) = fX 



— iiV^^ 



JL+^/s{2q/S)Yy_x+j/5{2q/S) - YL+^/s{2q/S)Jy_x+j/s{2q/S) 



JL+^/s{2q/5)Y^,s{2q/&) - YL+^/s{2q/S)J^/si2q/S) 



(14) 



To solve the dynamical part of dO), we seek a solution of the form: TZx,y{t) ~ T^x.yioo) — J2\ ''~y x^ '*^''i which leads 
to the following difference equation: 

^A , ^A , A , A , oA ,l+{y-x)5 

' y,2;+l y—l^x i^ y^x—1 i^ y-\-l,x 



2 (^A - {^^ ^}j r^y^, = 0, (15) 

with the boundary conditions: 

{L5 - Xq)rl^+L = and r^^ = 0, (16) 

Introducing the notation: E = ^\'^ , ( p^ admits as solution: r^y — AJy^x~E{2q/8) + BYy^x~E{2q/S), where A,B 
and the (relaxation-)spectrum {E} are determined from the boundary conditions ([l6|), which imply: 

" AJ-E{2q/S) + BY^E{2q/S) - ,^^. 

AJL-E{2q/S) + BYL-E{2q/S) = 0, ^ ' 

The only non-trivial solution of this system (for which A =/= and B ^ 0) requires 

JL-Ei2q/6)Y.Ei2q/6) - J^E{2q/5)YL-E{2q/5) = 0, (18) 



^ The case where Di — D2 = corresponds to the situation where the equations of motion of all the correlation functions 
close (see,(0), (M)) . For a translationally-invariant system with initial density p(0)and with Bi + B2 ^ Ci + C2, the density 

simply reads: (n,(t)> = a,+c,'-it+B,) + {pW " c+c.'-it+B,) ) 



or equivalently in terms of Lommel function 
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L-l.l-E 



{2q/S) = 0, 



(19) 



which admits L — 1 zeroes ||lJ,|l^,Q with degeneracy L. The latter are symmetricahy distributed around -j (which 
is an eigenvalue if L is even). To obtain the complete set of L(L — 1) + 1 eigenvalues, i.e., the complete relaxation- 
spectrum {Ei},i = 1, ... ,L of the string-function (which has not to be confused with the spectrum of the stochastic 
Hamiltonian H) , one has to take into account the eigenvalue qX = LS directly obtained from the boundary condition 
(Iq). Notice that in {Ei} the index i — 1, . . . , L labels the L distinct eigenvalues forming the relaxation-spectrum. In 
order to have some more insight into the relaxation- spectrum {Ei}, i = 1, . . . , L oi the string-function Sx,y{t) , we can 
take advantage from the fact that the following eigenvalue-problem: 



{E - n)Fn = V {Fn-i + Fn+i) ; (1 < n < L) 



(20) 



admits as eigenvalues the {L — 1) zeroes of the Lommel function: _Ri_i i_£;(2V^) — 119,0]. Therefore, choosing 
V ^ ^, the problem of determining the relaxation spectrum is reformulated as that of solving the eigenvalue-problem 
( pO| ) }A\J-)) — E\J^)), where A^ is a (L — 1) x (L — 1) symmetric (but not hermitian when q has an imaginary part) 
tridiagonal matrix and \!F)) is a (L — l)-components column- vector: \!F)) = {Fn=i F2 . . . F^^i)^ . The general form 
of the matrix Ai is the following: 
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q/5 3 q/5 

■•• ■■• ■■• 
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q/5 (i-2) 
q/5 



q/5 
{L-l)J 



(21) 



For small systems the [L — 1) distinct eigenvalues {Ei] of A^, can be computed analytically. For i = 6, with V = ^, 

we have {Ei{ — < 3, 3 ± J ^^'^^ ±V9+2av +4V l .^^j^gj-g .^^g gf^[\\ have to take into account the additional eigenvalue 

qX = L6. The spectrum depends on the size L of the system and this is in particular the case for E* = minE{E} = e^, 
the smallest eigenvalue, which governs the long-time behaviour of the system. For larger matrices we had to proceed 
numerically and for L 3> 1, cl — > £00, and E* is a constant: E* = eoo- For L = 6, we have the exact result 
eL=6 = 3 — w 5+4y +v9+24y +iv ^ xhis expression can be considered as an excellent approximation for systems of 
size L ^ 1 and in particular for eoo- 

Therefore, the long-time dynamics (of large systems, with L 3> 1) is governed by the eigenvalue {V = q/5): 



p* , ^o 5 + 4(1)^ + V9 + 24(1)^ +4(1)4 



(22) 



, I.e. 



qX* = E*5 + 7 = e^S + 7 



(23) 



This expression provides the inverse of the relaxation-time of the model under consideration. 

Having obtained the relaxation spectrum and the expression of r^ ^, the complete expression for the string- function 
follows as: 

S:,.,y{t) - 5.,,(oo) = ^^-^^B,e-(^'^+^)* [Jy^x^E,{2q/5)YL-EA'2q/S) - Yy^x^EX'2q/5)JL-E,{2q/5)] , (24) 

Ei 

Here for simplicity we consider the translationally-invariant (but not necessarily uncorrelated) situation, when 
Sx,y(t) = Sy-x(t). In this case, the coefficients AEi are obtained from the initial condition according to : 



L 

^E. = E K"'].j {Jn-Epql5)YL-E,{'2ql5) - Y„_EA^l/S)JL^EA2l/^)y (^nW " ^„(oo))/i-", (25) 



where A/" is an hermitian L x L matrix (see (pC 

To clarify this point let us introduce the following vectors of the Hilbert space C^ (with the usual scalar product): 

\S)) EE {{Sm - ^i(oo))a.-i . . . (^l(O) - 5l(oo))m-^))'^ (26) 

and 

/ {Ji.E, [2q/5)YL-E, {2q/5) - Yi-e, [2q/5)JL-E, {2q/5)) \ 

\Ve-)) = '■ (27) 

{Jl-i-e, {2ql5)YL-E, i2q/S) - Yl-i-e, {2q/S)JL-E, {2q/S)) 

V J 

In a vectorial formulation, the coefficients Ae are obtained from the initial state of the system according to: 

L 

\S))^Y.M\Ve,)) (28) 

Solving this equation, we formally obtain the expression for the coefficients Aej ■ 

L 

where Af is an hermitian L x L matrix which entries read: 

L 

[N],^^ = {{VeAVe^)) = E (J„-£.(2g/<5)lL-£,(2g/5) - y„_i=;.(2g/<5)Ji_£,(2(7A))* 

n=l 

X ( J„_B^ {2q/S)YL-E, {2q/S) - Y^-e, {2q/6)JL-E, {2q/6)) (30) 

With help of the expression of the string- function (p4h, we can compute the exact density of particles at site x: 



{nx{t)) = (31) 

In the non-translationally-invariant situation, we would proceed in a similar manner, but we would have to work 
with vectors of the Hilbert space C'^'-^"^-'"'"^. In this case we would have to take into account the degeneracy of 
the eigenvalues of M in order to compute the i(L — 1) + 1 components Ae appearing in (pih. This is achieved by 
replacin g E j with Ei^d in (|2J), where the index d labels the degeneracy of the eigenvalues Ei, i = 1, . . . ,L. 

From (|3l|), we can also obtain the expression of the non- instantaneous two-point correlation functions {nx{t)nxo{Q))- 
We should take into account the initial state ( pi| ) |-P'(0)) = n^„|P(0)) instead of \P{0)). 

With (|24|), we can also compute the instantaneous nearest- neighbour (two-point) correlation functions: 

/ \^,N « + Sx.x+2{t) - a{Sx,x+lit) + Sx+l,x+2it)) .„„x 

{nxnx+i)[t) = , (32) 

Although the present approach could be formally extended to obtain the exact and closed equation of motion of other 
string-like functions, which (eventual) resolution would provide the exact expressions of all instantaneous two-point 
correlation functions; in practice, such equations turn out to be extremely difficult to solve. The only cases where 
the whole hierarchy of equation has been completely solved are the free-fermion models. On the discrete (and finite) 
lattice these solutions were obtained by Hinrichsen et al. [n3| and in the continuum limit (for the diffusion-coagulation 
AA — > A model) by ben- Avraham [|2| . Unfortunately it is known that the latter approaches cannot be extended to 
systems which cannot be mapped onto free-fermion systems |l3,n2| . Here we prefer to take advantage of the quantities 
(iT-xit)) and (nxUx+iit)) and Sx^y(t), which we can compute exactly to obtain approximative instantaneous two-point 
correlation functions of the systems obeying the constraints (0), which associated string-function Sx,y{t) obeys the 
equation of motion (m. 



For technical convenience, we consider the translationally-invariant situation (thus {nx{t)) — p{t)) and expanding 
the string function we have: 






y-X>Ji>j2>--->Jn 

From equation (pS) , it is possible to obtain exact expressions relating a two — point correlation function with known 
quantities and higher order correlation functions. As an illustration, let us first consider the case where y — x = 3. 
Equation ( p3| ) implies: 

{nxnx+2){t) {nxnx+in,j:+2){t) = -^j4^ - -7- H 7 2{nj;nx+i){t) (34) 

For J/ - a; = 4, (||) with help of (||) leads to: 

/ \U\ f^\ I \(+\ Sx^:c+4{t)-2aS^j:^a^+3{t) /ON, 2 2ap{t) 

{na:nx+3)(t) ^ I ~ ) {nxnx+inx+2nx+3)[t) = j-^ h ( - 1 h {nxUx+ijit) (35) 

This procedure can naturally be continued for every two-point correlation functions. Therefore, for two-point corre- 
lation functions {nxnx+r){t) of sites separated by a distance r, using recursively the relations previously derived for 
{nxnx+r-i){t), {nxnx+r-2){t), ... we obtain an equality relating {nxnx+r){t) and an unknown r -I- 1-point correlation 
function {{uxUx+i ■ ■ ■ nx+r-inx+r){t)) to a combination of known quantities, as in ( |34| ) and ( |35|) . 

It is therefore possible to obtain approximative expressions for the correlation function within the truncation 
approximation (for r even): {uxUx+i ■ ■ ■nx+r^inx+r){t) ~ {nxnx+i){t) . . . {nx+r-inx+r)it) = [("-x^x+iXO]'"^^ and 
{uxUx+i ■ ■ ■nx+r~^inx+r){t) c^ [(n^rn-x+i) (t)] ''"~^)/^/9(i) , for r odd. 

Within this mean-field-like approach, we obtain the following approximative expressions for the two-point correlation 
functions: 

, w.x Sx,x+3(t) (^ , bp{t) \ 3ap{t) /aN2 
{nxnx+2)(t) :^ ^ I 2 + — — 1 {nxnx+i}(t) (^-j 

, w,^ Sx,x+4{t) -2aSx,x+3{t) I (b\^ , \t4-\\ I \u\ '^.apit) /a\2 
{nxnx+3)[t) ~ -— + II- I -I {nxnx+i){t) I {nxnx+i){t) — + (^-j 

; (36) 

To conclude this section, let us comment on this recursive procedure. 

First of all, the recursive character of the method appears through the repeated use of (33) and of the relations 
obtained for the other two-point correlation functions. The advantages of this recursive mean- field-like method with 
respect to the traditional mean-field are the following: 

(i) The procedure of truncation appears at the level of the three-point correlation functions. 

(ii) This approach is based on the explicit knowledge of the quantities Sy-xit), p{t) and {nxnx+i){t). 

(iii) This method does not give rise to non-linear partial differential equations and/or to non-linear self-consistent 
equations which are generally difficult to solve and which appear from traditional mean-field methods. Conversely, 
the approach presented here gives directly (after the truncation procedure) access to the (approximative) expressions 
of the two-point correlation functions. 

It follows from the exact expression (p4) of Sx,y{t) that, for the models under consideration in this work, the related 
string-function approaches its steady-state exponentially fast, with an inverse relaxation-time given by (|23|). This 
result is valid for an arbitrary initial state: the effect of the initial condition only appears through multiplicative 
factors Ae- In the translationally-invariant situation we have the coefficients (p5|); other initial conditions do not 
affect the exponential nature of the relaxation ( |24| ) with the inverse of relaxation-time (|23|). 



V. A MODEL WHICH CANNOT BE SOLVED DIRECTLY FROM THE CONVENTIONAL IPDF METHOD 

In this section we consider a model which cannot be solved directly by the conventional IPDF method. A brief 
account of the study of this model has recently been reported in pl| ] . Here we complete and develop this preliminary 
work. 

The dynamics of the model under consideration takes place on a finite and periodic lattice. When a particle and 
a vacancy are adjacent to each other, a branching reaction can take place and the particle A can give birth to an 
offspring {A% -^ AA and 9 A -^ AA) with rate FjJ = Fq}; another possible reaction is the death of the particle 
(^0 —>■ 00 and (dA -^ 00) with rate F?[5 = Fj]?. When two particles are adjacent, they can coagulate {AA -^ A9 
and AA — > 0.4) with rate F}° = F'j'}. In addition, when two vacancies are adjacent, a particle can appear {birth 
process, 00 -^ A9 and 00 -^ 0.4) with rate FJq = Fqq. The dynamics of this branching- coagulation with birth and 
death processes (BCBD) model can be encoded by the following reactions: 

00 -> (bA and 00 -^ A$ with rate: F°J = FJ° 
A9 -^ AA and 9A -^ AA with rate: F},^ = FJ} 
AA-* A^ and AA^ ^A with rate: FJ? = F?} 
AID -^ 00 and IDA -^ 00 with rate: F?° = F™ (37) 

It should be emphasized that in this model, the effective motion of the particles is realized by successive processes of 
branching, coagulation, birth and death on neighbouring pairs of lattice sites, without explicit diffusion process. 

The system described above can be viewed as an epidemic model where particles can spontaneously ap- 
pear/disappear, have an offspring and coagulate. It can also be viewed as a generalization of the voter model [0, 
where the presence/absence of particle is associated to an opinion (yes/no) and each site is associated to an human 
being. According to the dynamics of the model, each individual changes his opinion at a rate proportional to the 
opinion of his neighbours. 

For the model under consideration Ai = As = A = Fjj], Bi = B2 = S = F}J - Fjg, Ci = C2 = C = F^j] + Fjj], 
Di = D2 = D = Fj'q + FJq — (F}J + Til). If £* = 0, for the translationally-invariant situation with an initial density 

p(0) of particles, we have (when S 7^ C, see footnote (2)) {n^{t)) = ^^ + (p{0) - ^^b) e-^^^'"^)*. 

In this section we solve, with some restrictions on the reaction-rates, the above-mentioned model when D ^ 0, i.e. 
in the case where the equation of motion of the correlation functions of the model give rise to an open hierarchy (Q). 
It has to be stressed that this model can be casted into a free-fermion form only when FJq = F^q and F}^ = Too (see 
PI for a complete classification of free-fermion systems) . Furthermore, this model cannot be solved (directly) by the 
traditional IPDF method (not applicable |l^-|lj] in the presence of the processes AID -^ 00 ; IDA -^ 00 and in the 
absence of processes .40 -^ IDA ; 0.4 -^ AID; the latter should occur with the same rate as the coagulation rates 
|jll|-[l4|). The idea is to choose suitable a and b to close the equation of evolution of Sx,y{t). This is achieved by 
solving m) and thus imposing the following condition: 



F, 



= 1 + ^ > 1 (38) 



and the reaction rates fulfill: 



F}? = F?} > 0; 2FJS = 2^1 > F}J = T^ > 0; and F?S = F™ = ^"^'^"g^ ' ''^ > (39) 

We will see that the case treated in this section (with the constraints (|3t^)) can be obtained from the reversible 
model of diffusion-coagulation with input of particles (RDCI model), solved in the next section, via a similarity 
transformation. In fact, in section VII, we investigate for a local similarity transformation which would map the 
general Sx,y{t) function onto the empty-interval function (with a = b) and the RDCI model onto the present BCBD 
model. In the sequel we show that such a mapping exists and establish that the approaches followed in this section 
and in section VII are equivalent. 

For the model ( p7| ) with the restriction (^) and from the definitions (|o|), we have a = pi " (Fqq — F}q), /3 = 

1/2 
— |(Fqq — F}J) and thus a(3 < 0, fii ^ fi2 — f- ^ —i{sgna) ^ and q = ija/?!^". We also have 6 — 2bA/a > and 

-f^2{B + C)-d (Note that because of {^, < \q\/S < 1/2). 

The subcase F};^ = Fj^ implies a = (3 = B = D = and we recover (for C ^ 0) {nx{t)) = '^-s...+i{t) ^ 

u + (K(o))-u) 



Hereafter we focus on the more general situation where ( |39| ) are fulfilled with FjJ ^ Fqq, and thus a ^ 0, /3 7^ 0. 
The stationary expression of the string- function for this model is given by the expression (^^ . With help of the 
formula (^l|) and the ratio (g^) , we obtain the following expression for the stationary density of particles in the system: 



(na;(oo)) ^ -\a~ [i 



JL+-,,s{'iq/5)Y^+^,s{2q/5) - YL+^/si2q/S)J,+^/s{2q/S) 
JL+y/s{2q/S)Y^^s{2q/S) - YL+^/s{2q/6)J^/s{2q/6) 



(40) 



Similarly, with help of (^l|) and ( |38| ) , we obtain the stationary expression of the instantaneous nearest-neighbour 
correlation functions: 



(n^n^+i)(oo) 



2ab{nx{oo)) ~ a? 

62 



+ (^f 



JL+^,s{2q/5)Y^+^,s[2q/5) 



Yl 



+7 



,s{2q/5)J2+^,s{2q/5) 



JL+^/s{2q/S)Y^/s{2q/S) ~ Yl+^ /s{2q/ 5)J^/s{2q/ 5) 



(41) 



From (|4C|) and ( [4l|) one can check that the system under consideration is a correlated system of interacting particles. 
In fact, one can see that {nxnx+i){oo) ^ {nx{oo)) {nx+i{oo)) — {{nx{oo)))'^ , despite the fact that both steady-states 
( [40| ) and ( [4l|) are translationally-invariant, the stationary distribution is correlated which is due to the interacting 
character (hard-core) of the particles. We should emphasize that the presence of correlations in the stationary 
distribution is specific to the class of models considered here which cannot be mapped onto free-fermion systemsn (see 
also the model of the next section) . 

To study the dynamical properties of the model, we need the relaxation spectrum of the string-function. As 
established in the previous section, the latter is obtained as the set of zeroes of the following Lommel function: 
Rl-1 i-Ei2i\q\/5) = 0, where E = "^ ^^ . Solving the associate eigenvalue-problem (20) (in this case, the matrix A^, 
see (|2l[), is not hermitian). 

For small systems the (L — 1) distinct eigenvalues {Ei\ of (pl|) can be computed analytically. For L = 6, we have 



{E}= 3,3± 



5+4y^±V9+24y^+4\7T 
2 



For larger matrices we had to proceed numerically. Our analysis (based on 



the spectrum of large matrices, with L < 1000), shows that the spectrum {E} (and therefore {9A}) is always real and 
symmetric around L/2 which is an eigenvalue when L is even. The other eigenvalues are not generally integers, but for 
the central part of the spectrum (when eigenvalues which are close of L/2), the eigenvalues approach integer values. 
This is not the case at the extremities of the spectrum. In particular, the smallest eigenvalue E* = minE{E} is not 
an integer and depends on the size of the system: E* = e^ > 1. However, for L ^ 1, e^ ^ e^o, and E* is a constant: 

ForL: 



E 



> 1. 



6, we have the exact result eL=6 = 3--y/ ''+^'' +^"+^^^ +^'' , with 1 < e^^g < 3 

This expression can be considered as an excellent approximation for systems of size L ^ 1 and in particular for Coc 
Therefore, the long-time dynamics (of large systems, with L 3> 1) is governed by the eigenvalue E* — tL 



13. 



5+4y^+\/9+24y^+4yT 



I.e. 



q\* =E*S + -f^eLS + -f = 2 



^ 00^ 10 



^ lll^-L 00 ^ 107 



plO 
-^00 



(6^-i)(ri;j + ri^) 



> 2F}i > 



(42) 



The equation (E3) provides the inverse of the relaxation-time of the system. 

The dynamical approach towards the steady-state of the density is obtained from the dynamical expression 
the string- function, according to formula (^l[) and with the ratio (p8[), we obtain: 



{nx{t)) ~ {rixioo)) 



^^£.e-(^''+^^* [Ji-E,{2q/S)YL-EA'2q/S) ~ Y,^EA2q/S)JL-E.{'2q/S)] 



of 



(43) 



where the coefficients AEi have been computed for the translationally-invariant situation in 



d). 



^ In fact, for the free-fermion systems such as the diffusion-limited with pair annihilation and creation model |7| and the 
related diffusion-coagulation models [[Ll[-[L3|, it has been shown, for translationally-invariant systems, that {nxnx+r){oo) = 
({n,(oo)))2,Vr>0 



*As an illustration, for the case T^o = 3/10, TJi = 1/2, Til 



— 1/3, with the expression above, we obtain 



(analytically) eL=6 = 1.0823337683. For larger systems (L — 10,25,40,1000), we obtain numerically (with an accuracy of 



10" 



£10 — £25 — £40 = £1000 — 1.0823337697. 
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With help of (|^ and (^)_, we obtain the expression of the dynamical approach of the instantaneous nearest 
neighbour to its steady-state (f4l|): 



^^i;,e-(^'*+^)* [J2-E,i2q/S)YL-EA2q/6) - Y^-eA'^iI 5)Jl-eA'^<i/5)\ (44) 



Ei 



From (|24|), (|43|), ( |44|) and with help of (|33|)-(|3q), we can also obtain the approximative expression of the other two-point 
correlation functions. 

The result ( [43[) can be extended to the case of the non- instantaneous two-point correlation functions (n^ {t)nxa (0)) . It 

suffices to replace in ( [f3|) the coefficients Aei by those obtained from the initial condition ( {n?=x l*^ ~ ^"j (0))}fixo (0) ) ■, 

instead of Sy-x{Q)- 

Let us now mention the long-time behaviour of quantities computed above. For the sake of simplicity we consider 
the translationally- invariant situation, in the regime where E*t^ 1 {E* is the smallest eigenvalue (^2|)), we have: 

K(i)) - K(oo)) ^ (^) ^B.e-(^*^+^)* [Ji.E'{2q/S)YL-E'{2q/S) - Y,^E'{2q/S)JL-E'{2q/6)] (45) 






._, I). _(!)_, _(!),_ (I 



(46) 



It follows from the exact results (|4^)-(|46|), that the density {nx{t)) and the two-point correlation functions 
{nx{t)nxo{Q)) (non-instantaneous) and {nxnx+i){t) (instantaneous) approach the steady-state exponentially fast, 
as the string-function Sx.^jt), with an inverse of relaxation-time given by (P4). In addition, in the sense of the 
approximative scheme (|33|)-(|3q), these results are expected also to be valid for the more general correlation func- 
tions: (nxnx+r)(t) {r > 1) and for arbitrary initial conditions: the initial state only affects the multiplicative 
factors Ae of (^3|) and (Q) (for translationally-invariant systems the Ae are given by (P5|)). These statements 
are supported by the study of the subcase T\q = FJq (with C 7^ 0) which is solvable by conventional meth- 
ods and which dynamics is expected to be qualitatively the same as the more general case considered here (es- 
pecially for ITIq — rJol ^ 1). For this subcase the density has been computed previously (for arbitrary ini- 
tial condition) and the two-point correlation functions read: {nx{t)nx„{0)) = ^ + {{nxnxo){0) — ^) e~^'-^* and 

{n.ny)it) = {^f (1 - e-2^*) + {nxny){0)e-^^* + ^ ((^,(0)) + (n,(0)) - ^) (e-^^* - e-^^) , (y ^ x). 

VI. SOLUTION OF A REVERSIBLE DIFFUSION-COAGULATION MODEL WITH INPUT OF 

PARTICLES 

In this section we study a model of reversible diffusion- coagulation with input of particles which can be solved by 
the conventional IPDF method. Particles can jump (provided that the arrival site was previously empty) to the right 
and the left with rate F^J = FjJ > 0. We assume also that two adjacent particles can coagulate with the same rate 
^11 = r°} and that when a particle is adjacent to a vacancy, a branching process can occur with rate r}J = FjJ. In 
addition, when two adjacent sites are empty, a particle can spontaneously appear on a site (input) with (a finite) rate 
Fjg = FqJ > 0. The dynamics of this reversible diffusion- coagulation with input of particles (RDCI) model can be 
encoded by the following reactions: 

A9 ^ 9A with rate: r% = fJ? > 
00 -^ (bA and 00 -^ Ad) with rate: F[!J = FJ[J > 
A9 -^ AA and 9A -^ AA with rate: F}J = FJ} 
AA -^ A$ and AA -^ (bA with rate: F}? = F?} (47) 

This model can be solved by the conventional IPDF method: one can check that for the reactions (^) and ( [48| ) the 
constraints (ph are fulfilled with a = b. 

In order to apply the (conventional) IPDF method, we require that the coagulation and diffusion rate arc |0-n3: 

11 



^ 10 



pOl 

^ 11 



and 



plO 
^ 01 



r 



10 



(48) 



I) 



we have Fj'J 



pOl 

^ 11 



pOl 
J^ 10 



r|" > and are left with 4—1 = 3 independent reaction-rates for 



Therefore with 
the model (||). 

Before proceeding with the solution of this model some comments are in order. To our knowledge this model 
has been studied in the continuum limit on an infinite chain (for the translationally invariant situation) by Doering 
and ben-Avraham uM- The latter obtained, in this limit, the stationary concentration of particles, the stationary 
interparticle function and the relaxation spectrum as the zeroes of some Airy functions. In their work, Doering et al. 
considered an infinite chain with lattice spacing Ax (here Ax = 1 ) in the continuum limit {Ax — > 0) . On this infinite 
chain, the reactions occuring are the symmetric diffusion and coagulation with rates F^q 



plO 
^ 01 



pOl 

^ 11 



V 



plO 

^ 11 



r;ji = RAx. 



lattice of L 
> 0, 



2pio 

^J^ 01 



the symmetric branching processes with rate FjJ — Fq} = -^ and the input of particles with rate: FJq = 
In this section we want to solve the model (^) with the restrictions (^) on a finite and periodic 
sites. According to the definitions (10), for this section we have: ai — a2 = a = 2{Ci — ^i) = 
/3i = /32 = /3 = -2Di = 2Bi = 2(F?i + F}i - Fjlj), S^Ai+A2 = 2Fio > and 7 + (5 = 2(2Fi? + F}i). 

Hereafter we thus solve the model (E7h with the constraints (Es) , which is a model described by 4 — 1 = 3 independ 
paramaters (reaction-rates), for the case where the input of particles is non vanishing (i.e. S > 0) and D y^ 0. The 
case where (5 = have been extensively studied in the continuum limit lOJl^ as well as on discrete lattices ll^ : these 
models have been mapped onto free-fermion models p3|,p|,p[ . In addition, when D — 0, for the translationally- invariant 
situation and an initial density p(0) of particles, we have (when B j^ C, see footnote (2)) 



K(t)> 



A 



C-B 



P(0) 



A 



C-B 



,-2(C-B)t 



(49) 



D 



For the model under consideration here, we have /i = . /§ and q = \fa^. As a > and /? > 0, with /3 = 

i? = 0, we focus on the case where g ^ 0. One should be cautious with the fact that q and /i can be imaginary, in 
which case the sign ambiguity is fixed by requiring that qfj, — a. 

The equation of motion of the string-function associated to this model is thus of the form (|^) which has been solved 
in section IV. 

To compute the relevant physical quantities, we proceed as in the previous section. In fact the expressions of the 
density, correlation function and of the current of particles can be immediately obtained from (|40|)-(|4q) when ^ = 1. 
Hereafter for the sake of completness and clarity we quote these expressions: 



(n:r(oo)) 



Ai 



JL+^/s{^q/6)Y,+^/s{2q/6) - YL+^/s{2q/5)J^+^/s{2q/5) 
JL+^/s{2q/5)Y^/s{2q/5) - YL+^/s{2q/S)J^/s{2q/6) 



(50) 



{na:na:+i)ioo) = 2(n^(oo)) - 1 + ^^ 



JL+j/si'2q/S)Y2+^/s{2q/S) - ri+^/,(2g/5) J2+^/,(2g/5) 



JL+^/si2q/S)Y^^si'2q/5) ~ Yl+^ /s{2q/ 5)J^/s{2q/ 5) 



(51) 



We again note that the expressions dSQ-pTl) are independent of the site label x and therefore correspond to 
translationally-invariant stationary quantities. 

As for the BCBD model, from ( pO[ ) and ( ^l| ) one can check that the model RDCI is a correlated system of inter- 
acting particles, characterized by a translationally-invariant but correlated stationary distribution: {nxnx+i){oo) ^ 
{n^{(x>)){nx+i{oo)) = {{n^[oo))Y. 

To study the dynamical properties of the model, we need the relaxation spectrum {Ei\, i = 1, . . . , L of the string 
function. As established in the previous section, the latter is obtained as the set of zeroes of the following Lommel 
function: RL-i^i-E{2q/S) = 0, where E = ^-^, which are computed solving the associated eigenvalue-problem (|2C|). 
Notice that when q is real the matrix Ai (|2l| ) is hermitian, otherwise (g is imaginary) A^ is antihermitian. 

Again the spectrum {E} (and therefore {qX}) turns out to be real (even when M is antihermitian) and symmetrically 
distributed around L/2 which is an eigenvalue when L is even. 

When Fqq > FjJ -I- Fq} , then q is imaginary, with < -^ < ^ . When q is imaginary, the eigenvalues are not generally 
integers, but for the central part of the spectrum (when eigenvalues which are close of L/2), the eigenvalues approach 
integer values. This is not the case for the smallest eigenvalue E* = minE{E} which is not an integer and depends 
on the size of the system: E* = cl > 1- However, for L ^ 1, e^ — > Coo, and E* is a constant: _E* = too > 1- The 
expression (E2|) can again be considered as an excellent approximation for systems of size L 3> 1 and in particular for 
£00- We can then show that for q imaginary, the quantity el > 1. In fact one can compute (approximative) bounds 

for e^, with L ^ 1. This is achieved with help of (p2|), namely: 1 < eL=6 < 3 — a /2 + i 



13 . 



12 



< r, 



01 + rj}, q is real. The eigenvalues are still symetrically distributed around L/2 (which is still 
an eigenvalue when L is even). In this case, the smallest eigenvalue E* = e^ can be negative (e.g., eL=6 < if 



When rjo 



> 1.29907.. 



However, 

^10 



^J- 01 



< 



CL 



< 1 



in fact, 

10 I 



< 6^=6 ^ 1 )5 a-iid thus, when q is real, we have: 
> 0. Again, for large systems {L ^ 1) the expression 



2piu 

qX* = eL<5 + 7 = 2 [2^1 + r\l +"(eL - l)Tl"„] > 2(rio + eLVli) 
of Ci is well approximated by the exact expression ( |22| ) of eL=6- 

When rjj] = rj; + Tl\, then (with {^) D = B ^ q ^ and we recover (^), with B = 0). 

In definitive, the long-time dynamics of large systems is governed by the eigenvalue E* > according to (g2|) for 
L ^ 1. This means that the inverse of the relaxation-time of the system reads: 



qX* ^E*S + ^ = eL5 + j^2 [2ri° + T^ + {el - l)rl°] > 



The dynamical part of the above quantities are obtained similarly from 



setting - 



(52) 
1 in these expressions: 



24t)) - (n,(cx3)) = A*5I-4B,e-(^'^+^)* [Ji-E,i2q/6)YL-EA2q/6) - Yi^EA2q/S)JL-EA2q/6)] 



(53) 



("-x"-x+i)(0 - ("■x"-x+i)(oo) = {{nx{t)) + {n^+i{t)) - 2{n^{oo))} 

+ fi^Yl ^E. e"(^'*+^)* [J2-E, {2q/S)YL-E, {2q/d) - Y2-E, {2q/S)JL-E, {2q/S)] , 



(54) 



where the coefficients Ae have been computed in the translationally-invariant situation in (25)-(p0|). 

From equations (p3),(M) and (M), we can compute for this RDCI model the approximative expression all two-point 
correlation functions according to the scheme (|33| )- (|36| ) . 

The long-time behaviour of the quantities (pM-(p3) are also obtained as explained in the previous section and in 
particular the long-time behaviour of the density is given by ( [45| ) , where the smallest eigenvalue E* is the quantity 
obtained in (p^). 

As in the case of the model BCBD, and for the same reasons, the density and the two-point correlation functions 
(instantaneous and non-instantaneous) relaxe exponentially fast to the steady-state with an inverse of relaxation-time 
given by (H). 

Another relevant quantity which one can compute for this model is the so-called interparticle function px^y{t) 
[0,12|. The latter give the probability, a time t, for a particle at site a; — 1 to have as a next neighbour particle 
at a distant site y — x > 0. To obtain Px,yit), we set a — b ~ 1 and the string-function Sx,y(t) reduces to the 
empty-interval function |0,n2|: Sx,y{t) — {YVi=xi^ ~ "■i))(0 which is associated to the probability of having a 



sequence of holes starting from the site x and of length y — x. As shown by Doering et al. ]11| , |12| , it is possible to 
relate the density of particles, the empty-interval function Sx,y and the interparticle function Px,y(t) : Sx,y-\-i{t) — 
2Sx,yit)-\-Sx,y-iit) = ((1 -nx)-.-{l -ny_2) {uy^i - (1 -ny_i)ny)) (i). Therefore, with p{t) = 7; I]j("i(*))> for the 
translationally-invariant situation (the inverse of the p{t) measures the average distance between adjacent particles), 
we have: 



Px,yit) ^Py-x{t) 



[Sx.y + lif) — 2Sx,y{t) + Sx,y-l{t)] 



In particular in the stationary case, we have: 

Sx,y+l{00) - 2Sx,y{<X>) + Sx,y- 



Px,y{oo) = 



Jj. 



Pit) 



l(oo) 



(55) 



p(oo) 



.(f) 



pY, 



y—x+l+ 



2q 



Ci) + i^-% 



y-x~l+j 



l2q\ 
^ <5 ' 



2Y 



V-x-t-i\ S 



2£-j 



J, 



L+- 



■ ifWiCi) 



M^l- 



:(f)] 



Yr 



■(f)[^^(f) 



Pjl+ii 



2q^ 



.y-x 



YmC-§) 



HJy 



-x+i+'iiT) 



Ai ^Jy- 



-l+j(f) 



2.L 



-x+-iCj) 



Jl 



-iifWiCi) 



Mn+^(f)] 



y^+j(f)[j.(f) 



M+^(f)] 



(56) 



To conclude this section, let us point out the fact that the results (|5(]|)-(|5q) can be immediately generalized to 
systems which, in addition to the processes (H^), also include the (adjacent) pair-creation reaction: 00 — > AA, with 
rate TjJ. In fact, it suffices to replace Tjj] with Tjg -\- TjJ in (p0|)-(|56|) 



We also would like to emphasize the fact that the expressions (pC)-(Dg) are different from those obtained by Doering 
and ben-Avraham lllll who considered the continuum limit of this model on an infinite chain. 
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VII. SOLUTION, VIA SIMILARITY TRANSFORMATIONS, OF MODELS WHICH CANNOT BE 

MAPPED ONTO FREE-FERMION SYSTEMS 

In sections V and VI we have solved two different reaction-diffusion, (the BCBD and the RDCI), models which 
cannot be mapped onto free-fermion systems. It is therefore natural to wonder whether or not there exists a mapping 
between these models, transforming the empty-interval function (YVj^xi^ ~ "-j))(^) ii^to another string-function of 
the form (n?=x(l ^ t"j))(^)- ^^ ^^^^ section we study the class of models which can be obtained from the model 
of reversible diffusion-coagulation with input of particles (RDCI), analyzed in section VI through a class of local 
similarity trans form,ations. In so doing we will answer the following questions: 

(i) Does a similarity transformation exist which maps the RDCI onto the BCBD model and the empty-interval 
function (with a = b) onto a string-function (with - ^ 1) ? 

(ii) If so, does the mapping provide the solution of the model BCBD for the same constraints ( p9| ) considered in 
section V ? _ 

Let us consider the original stochastic Hamiltonian H, through the (local) similarity transformation B, define H as 

H = B-^HB (57) 

Because^ of the requirement that {d{t)){H,\P{0))) = (0(t))(7f, |P(0))), which implies that {x\Oe-"^\P({))) = 
{x\de-"*\P{Q)) = {x\dB-^e-"^B\P{Q)) = {x\OBB'^e-"^BB-^\P{Q)), it is clear^that under this transformation, the 
observable O, and the initial state |P(0)) transform according to: O = OB and \P{0)) = S^^|P(0)), where we assume 

homogeneous (uncorrelated, yet random) initial states with density < p(0) < 1 of particles: |P(0)) = I ('n^ 

In this section we focus on local transformations of the form P Jl^ , p^ 

L 

B = (g)b„ (58) 

i=i 

where bj denotes a 2 x 2 matrix acting at site j such that the stochasticity condition {x\^H = is fulfilled. In 
addition, in order to transform the empty-interval function into a more general string-function we want to consider 
transformations which map the operator 1 — Uj onto 1 — n-, = 1 — ruj, where r = -. Peschel et al. [H have shown 
that the transformation (p8|), with 

' 1 l-r 

3 



b, ^ ( , ) (59) 



satisfy this property. Through this transformation, according to ( p7| ) , the stochastic Hamiltonian related to the RDCI 

/ pOO pOO pOO pOO \ 

/ i,oo _oi ^10 ^11 \ 
model transforms into the following stochastic Hamiltonian —Hjj^i = 



pOl 

^ 00 


pOl 


pOl 

J- 10 


pOl 

i 11 


plO 
^00 


plO 

J^oi 


plO 
J- 10 


plO 

^ 11 


\ -"-oo 


pll 

^01 


pll 

^ 10 


pll / 
^ 11 / 



where the non-diagonal 



entries read: 



Tjo-nVr^n^^/r ; rJJ^O ; ^^l ^T^^> ^ {2nyTll){r - l)/r 

ni = rjo = mi - rlo)ir - 1) + rrj?]/r ; TJ} = T^ = VH/r ; T?? = 2(r - l)[(r - l)(r}i - rJJ]) -|- rFj^l/r 

r?} - r}? = [Tl^Ur If + (2 - r)((r - l)rij + rrj?)]/r (60) 

The uncorrelated and homogeneous, but random, initial state becomes: 

\pm = i\of j (61) 

One has to ensure that all the reaction rates appearing in ( |60|) are > 0, which require that r > and that < ^-^ < 1. 
Therefore we have the necessary condition: r > p(0) > 0. 
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According to the transformation (p8|), (159) the empty-interval function (11^=1: (l^"'j))(0 i^ mapped onto the string- 
function: {YljZ.^ (1 — nj)){t) = (nj=2;(l ~ '''^i))(^) ^iid therefore, the statistical quantities for the model described by 
the stochastic Hamiltonian H are obtained from the related quantities computed in section VI for the model RDCI 
as 

KW)5,i?(o)) = -KW>-ff:ip(o)> ; K«x+y)5|p(o)^(i) = ^Kri:E+y)fl-,iP(o)>(0 ;(y>o) (62) 



where, {nxlt))H^\p{o)) and {rixrix+i) H,\P(o)){t) have been computed in ( p3[ ) and (|54|). 
We will now consider two specific models described by (|60|). 
(a) To answer the questions (i) and (ii), we seek, through the mapping (pq) , (p9|) , a model of the BCBD-type 



and thus require, f^o = ^J? = Fj? = 0, as for the BCBD model considered in section V and infer from (|60|) 



{T\l - Tll){r - 1) + rVll = (r ^ 1), which implies 



pll _ plO 

tM ioo > (RX) 

11 I rio rio - ^ (.""^^ 



pll I T-io _ p 
-■- 10 + -•- 01 -*- 00 



Replacing the expression ( |63|) in (|60|), we obtain the reaction-rates of the new (BCBD) model in terms of the rates 
of the original (RDCI) model. In order to have a physical BCBD model, we have to require the reaction rates F's to 
be positive. We now take advantage from the fact that a version of the BCBD model has been solved in section V, 

where F™ — '^,{'^^w ~ ^\h)- I^ i^ therefore possible to check, from (pG) and (pSh, that this relation still holds in this 

case and parametrize the reaction rates ( pOj ) as follows: 

_ _ _ _ _ rio _ _ rio 

^00 ^00 

The requirement of positivity of these rates ( pif ) leads to 2T}^ > VYl- Thus, the reaction-rates ( p4|) describing the 

BCBD model from the RDCI model, are identical to those considered in section V for solving the BCBD directly from 

the generalized string function. It is therefore easy to obtain the density and the correlation functions of the BCBD 

model from the RDCI model inverting the relation (pi) and using (|6lj), (|62|). As an illustration, we consider the BCBD 

model with rates FJq = 1 ,F}J = 1/2 ,TYi = 2 and F^g = 3, which imply that r = 3. Such a model is thus mapped 

onto the RDCI model with the rates VH = 3 ,F}J = 3/2 ,Fj? = Fj? = 1. In this case the density of the BCBD 

model is related to the density of the RDCI model according to (|62|): {n.j.(t))^^^ ~ ^ ^ ^ — 

J 6 vtu; \ ^W(rjo=i x;j=i/2 ,rJo=2,r™=3;|P(o))) 

\{nx{t))f-^o^lj^ r"=3/'> ri"=ri"=iiP(o)^')" ^^ particular we have seen that the stationary density is independent 

of the initial state and thus, in this case, we have (see (E^) and (|50|)): {nx{(X)))^f^ ~^^ -^^ ~^^ = 

(1^00 = 1 '1^10 = 1/2 ■^11=2 ,rj(j=3) 
,1 (721/3)1-7 (V2i/3)-Y..i{V2i/3)J7. (72^/3) " 



3\'^^y°")'{rl°=3 ,r;j=3/2 ,rjt'=r;t'=i) — 3+3 



_fi_ 



^ 1 {V2i/3)Yi (■/2i/3)~Y^^ 1 {V2i/3)Ji (V2j/3) 



0.2401 (i> 1). 



We are now in a position to answer the questions (i) and (ii): 

We have shown that there exists a similarity transformation ^Sdj), ( pij ) that transforms the empty-interval function 
onto the generalized string-function and that maps the RDCI onto the BCBD model, with the same constraints of 
solvability that the constraints ( pilJ imposed in section V. We conclude that the present approach and the method 
devised in section V are equivalent. 

One additional comment on this equivalence is however useful at this point. Although both mentioned methods are 
equivalent, the method devised in section V is in a sense more convenient because it is direct: solving the equation of 
motion of the adequate (generalized) string-function, solves directly the BCBD model. Conversely, via the similarity 
transformation, we first solve the RDCI model, which is a task of the same difficulty as that of solving the BCBD 
model, and then find an adequate and non-trivial similarity transformation (where the new reaction-rates should be 
interpreted correctly in term of the original ones) . 

(b) Let us now consider a model which can be solved from the solution of the RDCI model via the similarity 
transformation (|5q) , (p9|) . The model under consideration is a reversible diffusion-coagulation with particles input and 
pair annihilatioiiplDCIPA) , which dynamics can be symbolized by the reactions: 

Afb <~^ fbA with rate: f% = fj? > 
00 -^ (bA and 00 -^ A$ with rate: r°J = rj[; > 
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A^- 


^ AA and %A 


-^ AA with rate: r\l = 


_fii 


AA- 


-> Afb and AA 


-> 0yl with rate: F}? = 


poi 
= J^ 11 


AA- 


-> 00 with rate: 


^00 

^ 11 





(65) 

This model(RDCIPA) can be obtained from the RDCI model via the similarity transformation (p8|), (p9|). Imposing 
rJJ = SFqq in ( |60| ) we get the following reaction-rates: 

flO_f01_plO/ f 11 _ f 11 _ 9^10 fOO — 9('7- - n fio 
^ 00 — ^ 00 ^ -'^ 00/ 'i -L 10 — J^ 01 — ^-L 00! J^ 11 — ^l' -l^M 01 

r}? = r?l = (2 - r)rj? + (r - i)rj°, rj? = rj? = ^r}i + r?J (66) 

For this (RDCIPA) model, we have three (positive) independent parameters: r > 1, Fqq > and Fj^ > 0. The 
positivity of the reaction-rates (|65|) and the physical meaning of the initial state requires the following constraints: 

r?i = (2 - r)f j; + (r - l)fJ!] > 0, r>l, < p(0) < r (67) 



Thus, for the model RDCIPA ( |65| ) described by the reaction-rates (p6f) with the restrictions (|67|), the density, the 
correlation functions can be computed from the results of the model RDCI according to (|62|), for homogeneous (but 
random) initial states described by (IgT 



VIII. PROPAGATION OF A WAVE-FRONT AND THE FISHER WAVES 

At the end of section II, we have stated that some reaction-diffusion models, are described at the mean-field level 
and in the continuum limit by non- linear partial differential equation of Fisher type (H) [pj . In this section we show 
that for some choices of the parameters (reaction-rates) the mean-field formulation of the models BCBD and RDCI 
gives rise to Fisher-type equations. Then, from the results obtained in section V and VI, we study the propagation of 
the wave-front from a microscopic point of view (in so doing, the correlation between particles are taken into account 
exactly). We show that the scenario predicted by Fisher's theory fails in one spatial dimension for the models under 
consideration. In all this section, we adopt the same notation as that introduced at the end of the section II. 

(i) For the model BCBD, setting 



_ (2Fi° + TTo Pii) + v/(Fii)2 + (F0°)2 + 4Fi°Fi° - 2FiiF°° 

'>BCBD = 



2[(r?S + rig)-(r}J + Fi?)] 



(68) 



and with the additional definitions: 



j^BCBD ^ 2^(rii)2 + (F;SJ)2 + 4Fj«Fi; - 2F}jF?o > , and k^^^ ^ 2 [f}J + F}° - (F™ + FJSJ)] > 0, (69) 



where the reaction-rates appearing in (|68|)-(|6S|) are those defined in 
(ii) For the model RDCI, setting 



_ 2Fi[j~Fli-HV(F}^)^+4F^iJF^y 

<Prdci — o rrio /'roi i rii\i 

^ l-*^ 00 ~ l-*^ 10 + -'^ io;J 

and with the additional definitions: 



(70) 



j^RDCi ^ 2^(r}i)2 + 4F10F1; > and k^^^' ^ 2 [(FJ? + F}J) - FJ"] > 0. (71) 

The reaction-rates appearing in ([70|)-([7l|) are those introduced in section VI. 

Under the conditions (i) and (ii), at the continuum mean-field level, we have for the model BCBD and RDCI, with 
l^Mpi^^ t) = p^j^p{x, t) — (j)i, where I = BCBD, RDCI, equations of motion which are Fisher's equations: 

d fc' 9^ _ _ 

-g-fiiF(.x,t) = -j--^AiF{x^i) + k[p'-j,.jp{x,t) - k^2ipiiFix,t)f (72) 

Assuming L to be even and relabelling the sites of the chain according to the shift: x — > x— ^, we consider an initial 
inhomogeneous configuration with 
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K(0)>-|^^^??^^'^^''^t~^/^'"^ (73) 

^ ^ '' 10, otherwise ^ ' 

We want now to compare the prediction of the mean-field theory with the results obtained directly from the microscopic 
results derived in section V and VI and thus compute the time-dependent position X{t) of the wave-front and its 
time-dependent width w{t). This is done according to the formulae [ [L7[ : 



With help of (|4^) and (|52|), and denoting with qi,Si,^i,Ei the quantities related to the model I e {BCBD,RDCI) 
and defined (computed) in sections V and VI, we obtain in the long-time regime {EiSi -t- 7;)i ^ 1 where the time 
scales as i ex L^: 



Xi{t) = ^2ui{E*Si + -fi)t [l + 0{e-iVi^y^')^ (75) 

and 



wi{t)^Jui{E;Si+-/i)t l + 0(e-iV^^^) , (76) 



where we have introduced the parameter ui = 2(e'^ + — it ~ ^i^)- 

From these exact results, it appears that the location of the wave-front moves as Vt. Moreover, in contrast to the 
Fisher's mean-field theory, the width of the wave-front broadens as yt- These results which have also been observed, 
in the continuum limit, for the one-dimensional reversible diffusion-coagulation (without input) model [|l7| , confirms 
that in one spatial dimension the mean-field Fisher's picture fails. In fact Riordan et al. M^ have argued on the basis 
of numerical computations for the reversible diffusion-coagulation (without input) model that in higher dimension 
(> 4) the latter model is in agreement with Fisher's mean-field predictions: the width of the wave-front does not 
broaden. Very recently other authors who studied the same model as Riordan et al. (in dimensions d > 1) came to 
completely different conclusions |23]. 

Furthermore, for Pjyjp, Fisher's equation admits two (homogeneous) stationary-states, namely Pjv/_f(°°) = it which 

is linearly stable p6| and p^j^(oo) = which is linearly unstable flH]. This implies that at mean-field level, we 
would have for the stationary density, p\.^p{<x>) = p5\j^(oo) -f 0; which corresponds to two possible steady-states: 

p\jp{oo) = TT + </"/, Pm_f(°^) ~ ^i- However, from the exact expressions of the stationary density (^0|) and (|5C|), 
we know that the models under consideration admit unique steady-states which do not coincide with the mean-field 
prediction. 

IX. SUMMARY AND CONCLUSION 

In this work we have extended the conventional interparticle distribution function (IPDF) method. We introduced 
a string-function which is a natural generalization of the empty-interval function employed in the IPDF method. 
We derived the (five) constraints for the equations of motion to close (see (||) and (||)). We solved the equation of 
motion of this string-function on a periodic and finite lattice for the general form of a class of models which cannot 
be mapped onto free-fermion systems and which so far (to our knowledge) have been poorly understood (see (24[)). 
Then we specifically studied two models: The first one, which is a model with branching, coagulation, birth and death 
processes (the BCBD model), can be viewed as a generalization of the voter model and/or as an epidemic model. 
The BCBD model is an example of model which cannot be solved directly by the traditional IPDF method. For 
this model, under certain restrictions on the reaction-rates (see (|39|)), the density, the non-instantaneous two-point 
as well as the exact nearest-neighbour (instantaneous) correlation functions have been analyzed: the steady-states 
(see (|40|)-(|4l|)) as well as the dynamical approach towards the latter has been computed exactly (see ([l3|)-( [44|) ). In 
particular the relaxational spectrum as well as the the inverse of the relaxation-time have been obtained (see p2|)). A 
similar analysis has been performed for a reversible diffusion-coagulation model with input of particles (RDCI model). 
The latter (with the usual restriction that the coagulation rate coincides with the diffusion one) can be studied with 
help of the traditional IPDF method (the string-function then reduces to the empty-interval function) . In addition to 
the above-mentioned quantities, which we were able to compute also for the RDCI model, we calculated the stationary 
interparticle- function (see (|5(|)). 
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On the basis of the exact results, we have developped an approximative recursive scheme that allows to compute 
the (other) instantaneous two-point correlation functions (see (|36|)). 

Studying these (BBCD and RDCI) models, we observed that the latter are characterized by a translationally- 
invariant stationary distributions, for which contrary to what happens to free-fermion systems correlations are present: 
(n^n^+i)(oo) ^ ((n^(oo)))2. 

Later we studied the solution of the RDCI model and its implications on other systems related via similarity 
transformations. In particular we considered a class of similarity transformations (see (pq),(p9h) which transforms 
the conventional empty-interval function into a more general string-function. In so doing we saw that it is possible 
to map the RDCI model onto the BCBD one , which turns out to be solvable (via the similarity transformation (pq), 
(p9|)) with the same constraints encountered in section V. We therefore conclude that the approaches of the section V 
and VII for solving the BCBD model are equivalent. However it has to be noticed that working with the generalized 
string-function as in section V gives naturally access to the solution of the BCBD model without requiring the solution 
of another (the RDCI) model. 

We also have identified a model of reversible diffusion-coagulation with particles input and pair annihilation (RD- 
CIPA) which can be mapped, for some choices of the reaction-rates (see (|66|)-(|67|)), onto the RDCI model. For this 
RDCIPA model all the quantities previously computed for the RDCI can be immediately obtained via the similarity 
transformation (see (|62|)). 

Finally we observed that on some parameter manifold, the mean-field approximation of the BCBD and RDCI 
models are described (in the continuum limit) by the so-called Fisher equations which predict that an inhomogcneous 
initial configuration will evolve without broadening of the wave- front in the density of particles. Computing the 
width of the wave-front which broadens as Vt (see ([76|)), we show that the Fisher's mean-field description fails in 
one-dimension. Another failure of the mean-field theory is observed when one compares the mean-field predictions 
for the steady-states of the density with the exact results which shows that these stationary states are in fact unique. 
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